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We investigate local three-body correlations for bosonic particles in three and one dimensions as a function 
of the interaction strength. The three-body correlation function g' ' is determined by measuring the three-body 
recombination rate in an ultracold gas of Cs atoms. In three dimensions, we measure the dependence of g^> 
on the gas parameter in a BEC, finding good agreement with the theoretical prediction accounting for beyond- 
mean-field effects. In one dimension, we observe a reduction of g' ' by several orders of magnitude upon 
increasing interactions from the weakly interacting BEC to the strongly interacting Tonks-Girardeau regime, in 
good agreement with predictions from the Lieb-Liniger model for all strengths of interaction. 
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Correlation functions reflect the non-classical nature of 
quantum many-body systems. They may be used to character- 
ize the latter when quantities such as temperature, density, di- 
mensionality, and particle statistics are varied in experiments. 
It is particularly instructive to monitor a system's correlation 
functions as the strength of particle interactions is tuned from 
weak to strong. A paradigm is given by an ensemble of bosons 
in one-dimensional (ID) geometry with contact interactions 
[1 j: For weak repulsive interactions, in the zero-temperature 
limit, the system is a quasicondensate with essentially flat 
particle correlation functions in position space to all orders. 
For strong repulsive interactions, the bosons avoid each other, 
leading to loss of coherence and strong increase of local 
correlations. In the context of ultracold atomic gases, with 
exquisite control over temperature, density, and dimension- 
ality 0, tuning of interactions is enabled by Feshbach reso- 
nances [3|. Local two- and three-body correlations in atomic 
many-body systems can be probed e.g. in measurements of 
photoassociation rates [4| and of three-body recombination 
processes [5, 6|, respectively. Non-local two-body correla- 
tions for atomic matter waves have been measured in atom 
counting ll7l- [T0l . noise-correlation ifTTHTJl . and in-situ imag- 
ing 1 14 1 experiments. Recently, also non-local three-body cor- 
relations have become accessible in experiments [T51 IT61 . 

Recombination processes are sensitive to the properties of 
the many-body wave function at short distances. In particu- 
lar, the process of three-body recombination, in which three 
particles collide inelastically to form a dimer, is directly con- 
nected to the local three-particle correlation function g' 3 ' = 
(\j/ if (x) i \j/(x) i ) /n 3 , which compares the probabilities of hav- 
ing three particles at the same position for a correlated and an 
uncorrelated system. Here, i// T and \(f are atomic field opera- 
tors and n is the density. The function g^' depends strongly on 
quantum statistics ll5l [T6l and temperature lfT7l[T8l . For exam- 
ple, in 3D geometry, statistics change the value of g( 3 > from 
zero for identical fermions to one for non-interacting classi- 
cal particles and to six for thermal (non-condensed) bosons. 
For non-interacting bosons statistical bunching is suppressed 
in a Bose-Einstein condensate (BEC), for which g^' = 1. In 



addition, interactions also have a pronounced effect on g^>: 
In a 3D BEC, quantum depletion due to quantum fluctuations 
reduces the condensate fraction by increasing the number of 
occupied single-particle modes. In this case, beyond-mean- 
field calculations [19] predict an increase of g' 3 ' proportional 
to the square root of the gas parameter (na^ D ) 1 ' 2 , where «3d 
is the 3D s-wave scattering length. This increase of g( 3 > has 
never been seen experimentally and is in stark contrast to the 
behavior of ID systems. In ID geometry, bosons with repul- 
sive interactions minimize their interaction energy by avoid- 
ing spatial overlap. For very strong repulsive interactions in 
the Tonks-Girardeau (TG) limit 1 1 , 20-23 1 a strong reduction 
of g^' with a y~ 6 scaling is predicted l24l . Here, 7 is the di- 
mensionless Lieb-Liniger parameter, which characterizes in- 
teractions in a homogeneous ID system 1 1 , 25 1. Recently, g( 3 ' 
has been calculated all the way from the weakly to the strongly 
interacting ID regime |26|. Experimentally, Laburthe Tolra et 
al. (6) have observed a reduction of g^> by a factor of about 
7 (5) for a weakly interacting gas of Rb atoms with 7 = 0.45. 

In this work we experimentally determine g@> in 3D and in 
ID geometry using a trapped ultracold gas of Cs atoms with 
tunable (repulsive) interactions. For a BEC in 3D geometry 
we find clear evidence for an increase of g' 3 ' with increasing 
interaction strength, in good agreement with the prediction of 
Ref. 1191 . In ID, for which we can tune 7 from zero to above 
100 1 23 1, we determine g^> in the crossover regime from 
weak (ID BEC regime) to strong interactions (TG regime). 
Here our data agrees well with the prediction of Ref. [26 1 . For 
strong interactions in the TG regime, our measurements show 
that g( 3 > is suppressed by at least three orders of magnitude. 
For high densities and strong interactions, we observe a rather 
sudden increase of three-body losses after long hold times in 
the trap. Understanding the behavior of g^ at short and long 
times is an important step towards understanding integrability 
and thermalization in ID systems Il27ll28l . 

A three-body loss process |[3][29) consists of the collision of 
three particles, the formation of a dimer, and the release of the 
dimer' s binding energy typically sufficient to allow both, the 
dimer and the remaining particle, to escape from the trap. The 
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FIG. 1: (color online) (a) Relative atom number N(t)/N(0) vs. hold 
time t in 3D geometry: BEC (squares and circles) and thermal gas 
(diamonds) for a^u = 101(2) uq, 386(3) uq, and 386(3) arj> respec- 
tively. The dashed lines are fits to the data based on the loss equa- 
tion (see text). The solid lines are linear fits that include the data 

(3) (3) 
from 100% to 85%. (b) The ratio of correlation functions g lh /^bec 

as a function of a^u (experimental data: circles; prediction 1191 : 

squares). All error bars reflect the lo" statistical uncertainty. 



loss, assuming negligible one- and two-body loss, is modeled 
by the rate equation h = —aK^ 3 'g^ 3 'n 3 . Here, we have explic- 
itly split the loss rate coefficient aK^'g^ 3 ' into its three con- 
tributions. The parameter a = 3 describes a situation where 
exactly three particles are lost in each recombination event. In 
principle, secondary losses 1 30 1 could modify its value. How- 
ever, in the following we will be interested in relative mea- 
surements of aK^g^ 3 \ which are only weakly dependent on 
the precise value of a [25 1, allowing us to neglect a possible 
deviation of a from the value of 3. The parameter K^> con- 
tains the effect of few -body physics on the loss process Q. It 
depends on the probability of dimer formation (a process that 
can be strongly enhanced near Efimov resonances [31 1) and 
generally varies strongly with «3d [29, 32-351. For 030 much 
larger than the range of the scattering potential, K^ 3 ' shows a 
generic a 4 D scaling. Contributions of many-body physics are 
contained in the three-particle distribution function g( 3 'n 3 . In 
what follows, we aim to measure g( 3 ' as a function of a^ both 
in 3D and ID geometry. 

We determine K^'g^ 3 ' from measurements of the decay of 
the total number of atoms N(t) in our trap (5j|35), which obeys 
the loss equation N = —3K@>g@> fn 3 (r)d 3 r. Figuresllla) and 
|2lc) show typical atom number measurements for 3D and ID 
geometry. The data in 3D geometry is well fit by solutions 
to the loss equation. The determination of K^ 3 'g^ 3 ' depends 
critically on an exact knowledge of the atomic density profile 
«(r). In particular, particle loss and loss-induced heating of 
the sample [35 1 can modify the density profile in a non-trivial 
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FIG. 2: (color online) (a) Sketch of the experimental setup: A 2D 
optical lattice traps atoms in an array of ID tubes, (b) Example of 
a computed atom number distribution Nu (see text), (c) The rela- 
tive atom number N(t)/N(0) as a function of time t in ID geome- 
try: squares and circles correspond to a^u = 23(1) ciq and 568(3) a$ 
with initial densities of 4.5 ,um _1 and 1.7 /irn"' at the center of the 
center tube, respectively. The solid lines are linear fits to the ini- 
tial slopes, (d) The relative atom number N(t)/N(0) in ID for fixed 
a^u = 568(3) ao and for various values of y as the ID density is 
changed: y = 12 (circles), y = 13 (triangles), and y = 14 (squares). 
The solid (dashed) lines are linear fits to the data points for short 
(large) times to guide the eye. 



way. Also, on long time scales evaporative losses might start 
to play a role. To avoid these complications we restrict our- 
selves to short time intervals, during which not more than 15% 
of the atoms are lost, and we determine the slope N(Q) from 
a linear fit to the data. We determine Jn 3 (r)d 3 r from a mea- 
surement of the total atom number N and the trap frequencies 
(0 X yz usm g interaction dependent models for «(r) |25l . We 
find that the linear approximation underestimates K^ 3 'g^ 3 > by 
approximately 12%, however, the data analysis is greatly sim- 
plified, especially in ID. Finally, a comparative measurement 
of K( 3 >g( 3 > allows us to eliminate K^ 3 \ as explained below, 
and to determine g( 3 > in 3D and ID geometry. 

Correlation function in 3D: We measure K^ 3 'g^ 3 ' for both 
a non-condensed thermal sample and a BEC as a function of 
a3D- For the thermal sample we start with typically 3.5 x 10 s 
Cs atoms at a temperature of T « 200 nK. The peak density 
is about n = 1 x 10 14 ctrr 3 . In the BEC [35] [36] we have 
about 9 x 10 4 Cs atoms without any detectable non-condensed 
fraction at about «o = 5 x 10 13 cm" 3 . We tune 030 in the 
range from 50 «o to 800 «o by means of a broad magnetic 
Feshbach resonance [35, 53 («0 is Bohr's radius). The mag- 
netic field gradient needed to levitate the atoms against grav- 
ity (35'l introduces a slight (less than 5 oq) variation of «3d 
across the samples. We determine N by means of absorption 
imaging after a variable hold time t and 50 ms of expansion 
in the presence of the levitation field. We note that we do not 



observe the appearance of any non-condensed fraction in all 
measurements using the BEC. Figure [Tib) displays the ratio 
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mal sample and the BEC as a function of a^. Here we have 
made the reasonable assumption that K^' is independent of 
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pected value of 6 for weak interactions 1 5 1 , but then exhibits a 
pronounced decrease as a-yo is increased. For comparison, we 
plot the prediction of Ref. Ifl9l 
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We note that the density no enters into this equation as a mea- 
sured quantity. In general, we find good agreement between 
the experimental and the theoretical result, establishing our 
measurement as a clear demonstration of beyond mean-field 
effects on g^> in 3D bosonic quantum gases. 

Correlation function in ID: Figure [2] (a) illustrates our ex- 
perimental setup to generate an array of ID systems. We load 
a BEC of typically 8 x 10 4 atoms within 400 ms into approx- 
imately 5000 vertically (z-direction) oriented tubes that are 
formed by two horizontally propagating, retro-reflected lat- 
tice laser beams. Each tube with index (i,j) in the x-y-plane 
has a transversal trapping frequency of co± = 2% x 12.2(5) 
kHz and an aspect ratio 03±/(0 Z of approximately 800. The 
transversal motion of the atoms in the tubes is effectively 
frozen out as kinetic and interaction energy are much smaller 
than h(Q±. We adjust aj£> in 100 ms to its final value. Af- 
ter time t we turn off the lattice potential and determine the 
total atom number N(t) by absorption imaging in a time-of- 

flight measurement. 
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we calculate 



In order to determine g\ D 
^3)„(3)>, _ „(3)/gg) t Here, it is not obvi- 
ous that few-body physics is not affected by the confinement 
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and that hence Zf 1D and K^ cancel each other. Neverthe- 
less, it is reasonable to assume that K™ is not significantly 
changed by the confinement as long as the confinement length 
ax = v h/(ma>j_) is larger than the extent of the dimer pro- 
duced in the recombination event and the range of the scat- 
tering process, which are both of order of a^. Here, m is 
the atom mass. We choose a moderately deep lattice potential 
with a± m 1500 ao and restrict 030 to 030 < 800 ao- In partic- 
ular, we avoid the confinement-induced resonance condition 

fl 3 D«fl± E2H23). 

(3) (3) 

The main difficulty in the determination of K lD g 1D comes 
from the fact that the initial atom number of the tubes varies 
across the lattice as a result of the harmonic confinement. We 
choose to always load the lattice in a regime of weak repul- 
sive interactions such that almost all ID samples are initially 
in the ID Thomas-Fermi (TF) regime ll39ll . The local chemical 
potentials ji, >j, the total atom number N, and the chemical po- 
tential fi are then unambiguously related, and we can directly 
calculate the initial occupation number N{< for each tube (i, j) 
fll25ll and Fig.[27b)). The variation in Ntj results in a consider- 
able variation in the type of density profile for each of the ID 



systems after the strength of interactions is increased to the 
desired value: Some tubes remain in the ID TF regime, while 
others are now in the TG regime. For tubes that are in the 
weakly interacting regime we determine the ID density «i D 
numerically by solving the ID Gross-Pitaevskii equation. For 
the TG regime the density profiles are determined following 
Ref. |39|. In general, we find good agreement when we com- 
pare the numerical results to integrated density distributions 
from in-situ absorption images. For the interaction parameter 
7 we take a mean value that is calculated as an average over 
all local jij at the center of each tube (?', j) weighted by Nu 

ma. 

(3) (3) 

As before we determine ^ 1D g 1D from the initial slope of 
the loss curve as shown in Fig.[2jc). In Fig. [3ta) we compare 

(3) (3) 

the data that we obtain in ID geometry to our data for K^ D g 3D 
for a 3D-BEC as we vary 030. We note that the BEC data is in 
good agreement with previous three-body loss data on thermal 
samples when one takes into account the combinatorial factor 
3 ! = 6 13T1I351 . In particular, the 3D data follows the universal 
scaling law TO 3 ) ^ a ^ f or sufficiently large a^o E9l [3214351 . 
We exclude data points affected by the presence of a narrow 
Feshbach resonance in the vicinity of 030 = 150ao |40|. Note 
that in the range from 030 ps 10ao to 030 sw 850ao three-body 
losses in 3D increase by nearly 3 orders of magnitude. This 
behavior is in stark contrast to the measurements in ID. In ID, 
we observe a reduction of K^g^ by approximately a factor 
of 2 upon increasing 030 over the same range of values. In 
fact, for A3D > 200flo our measurement only gives an upper 

(3) (3) 

bound on ^ 1D g 1D as losses become so small that we have 
difficulty in determining N(0). Note that tunneling between 
tubes (on a timescale of 1 s for the parameters of our lattice) 
sets an upper bound for the timescale for which the tubes can 
be considered to be independent and hence fully in the ID 
regime. 

In Fig. gb) we plot K®g® / (A§*§ ) « g$ as a function 
of 7 A striking decrease by 3 orders of magnitude from the 
value 1 at 7~ 0.03 to 10~ 3 at 7 sa 50 can be seen. We compare 
this result to the predictions based on the Lieb-Liniger model 
of interacting bosons in ID: In the weakly interacting Gross- 
Pitaevskii regime (7C 1) the Bogoliubov approach yields 
g( 3 )(y) ~ 1 —6^/y/n, while in the TG regime, 7» 1, g' 3 ) can 
be expressed through derivatives of the three-body correlation 
function of free fermions, giving g^< — 167T 6 /(15t 6 ) J24l . 
Cheianov et al. f26l have recently calculated numerically g^> 
for all strengths of interactions within the Lieb-Liniger model, 
providing an interpolation between the weakly and strongly 
interacting limits (red continuous line in Fig. [3jb)). We find 
very good agreement between the result of our experiment and 
the theory that is valid for all strengths of interactions. This is 
the central result of this work. 

Finally, for large values of «3d and no, and for long hold 
times in ID geometry, we find a surprisingly sudden increase 
of losses as shown in Figpld), accompanied by a rapid in- 
crease for the expansion energy in the longitudinal direction 
(data not shown). The onset of increased losses shifts to later 
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FIG. 3: (color online) (a) Three-body loss coefficient K^ >g( > vs. 
«3D for a quantum degenerate gas in 3D (squares) and in ID (cir- 
cles). The line gives the A'' 3 ' = Cha^/m scaling in the universal 
regime in 3D with C = 67.9 |29ll32lf35l . The error bars of tf< 3 V 3 ' 
reflect the 1 o~ statistical uncertainty of the linear fit. (b) The mea- 

(3) 
sured correlation function g^ vs. y in ID geometry (circles). The 

(3) (3) 

crosses indicate the values for g lr , corrected for the variations of g™ 

with «3d as in Eq. {TJ (see text). The dashed and dotted lines are ap- 
proximate analytical solutions for ycl and ys> 1 from Ref. 1241 . 
The solid line is the prediction from Ref. 1261 . 



times with decreased density in the tubes, i.e. increased y, 
and it is rather sensitive to the precise value of y. We believe 
that the ID tubes suffer from a recombination-heating induced 
breakdown of correlations: For sufficiently large values of «3d 
the binding energy of the weakly bound dimer produced in the 
recombination process becomes comparable to the trap depth 
(here h x 45 kHz). This leads to a positive feedback cycle in 
the many-body system in which three-body losses lead to an 
increase of temperature ll3~5ll and thus of g^> ifTTl . which in 
turn increases three-body losses. 

In summary, we have measured the local value g( 3 > for 
the three-particle correlation function for quantum degener- 
ate gases in 3D and ID. In 3D, increasing interactions deplete 
the condensate and increase the value of g^ in accordance 
with beyond mean-field calculations. In ID, we observe a 
strong suppression forg' 3 ' by 3 orders of magnitude as the TG 
regime is entered. The accompanying suppression of three- 
body losses is crucial to the study of strongly interacting mat- 
ter in and out of equilibrium in ID lr23ll27llJ8ll4T1 . 
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SUPPLEMENTARY MATERIAL 
Trap parameters 

In 3D geometry, we measure the atom loss in a crossed 
beam dipole trap with one horizontal and one vertical laser 
beam. The horizontal trap-frequencies (0 x . y and the verti- 
cal trap-frequency (0 Z vary for the different measurements. 
The data sets in Fig[TJa) are taken with trap frequencies 
C0 A , y , z = 2k x (29(1), 80(2), 74(1)) Hz for thermal atoms and 
with ffl w = 2% x (11.8(1), 17.9(3), 13.5(1)) Hz for a BEC. 
The data sets in Figpja) are taken at trap frequencies of 
C0 Xi y, z = 2nx (10.5(8), 17.4(1), 13.9(1)) Hz for a BEC. 

In ID geometry, we use a crossed dipole trap in addition to 
the 2D optical lattice potential to adjust the atom number dis- 
tribution over the tubes. We choose two settings with global 
trap frequencies 0) r . ViZ = 2% x (9.7(2), 11.4(2), 14.5(1)) Hz 
and 2% x (13.1(2), 17.7(2), 17.5(2)) Hz. 



Atom number distribution over the tubes 

We calculate the initial occupation number Ntj for tube 
(i,j) from the global chemical potential ji. For weak repul- 
sive interactions during the loading process almost all tubes 
are in the ID Thomas-Fermi (TF) regime with a local chemi- 
cal potential jU,j at the center of each tube 



IMj = M 



1 



>(A/2)>^ + </), 



where m is the atomic mass and X = 1064.5 nm is the wave- 
length of lattice light. We calculate ji from the condition 
N = LijNijin) with the N Lj given by |Q] 
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where g\o is the ID coupling parameter [2| 



g\u = 2fico±a^ 1 — 1.0326 
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Determination of y 

We determine the mean interaction parameter y from the 
local parameters Yij at the center of each tube (i,j) 



Yi, 



mgip 
% 2 n}f 



xL N ur,j- 



Here, n\j is the ID density at the center of the tube (i,j). 
Note that this gives a lower estimate for y. Averaging y over 
the density profile along each tube gives a slightly larger y by 
a factor 1.5 for a ID TF density profile and a factor 1.27 for a 
TG density profile. 
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FIG. 4: (color online) The quantity JnfJr)dr for an occupation 
number N(j = 15 as a function of the scattering length 030 for the 
given trap parameters of our experiment. The various curves cor- 
respond to the different approximations: gaussian solution (green 
dashed line), TF solution (blue dotted line), TG solution (black dash- 
dotted line), numerically solved GP-equation result (black stars), and 
Lieb-Liniger solution with local density approximation (blue cir- 
cles). For the data analysis the continuous red line is used. 



Density profiles 

The density profiles for the individual tubes with index (/, j) 
depend strongly on the strength of interactions and the occu- 
pation number Nfj. Fig. [4] compares the results for the inte- 
grated density profiles jnfAr)cPr using the different approx- 
imations to calculate the profile (gaussian, TF, TG, numer- 
ically solved GP-equation, and Lieb-Liniger solution within 
the local density approximation |3|) for the specific case of 
Njj = 15. For our analysis of the experimental data we use 
the GP result for weak interactions and the TG result [4] for 
strong interactions (continuous red line). 



Secondary loss processes 

Here we estimate the deviation Aa from a = 3 in the rate 
equation h = —aK^g^'n 3 due to secondary loss processes 
15] ■ Within a simple simulation, we determine an upper 
bound for the correction to the data of the 3D loss experiment 
of Fig. f]7b) and show that secondary loss processes cannot 
explain our results for g th /g B ^ c - In fact for our experimental 
trap parameters and atom numbers secondary processes would 

(3) (3) 

result in an increase of g th /g B ^ c with increasing interaction 
strength, in contradiction with the observed behavior. 

A secondary collision is caused by the collision of the dimer 
and/or free atom from a three-body recombination process 
with other atoms while leaving the trap, triggering additional 
losses. We estimate the average number of secondary colli- 
sions in our experiment by determining numerically the colli- 
sional opacity («/) a for the products of a three-body recom- 
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turn of the dimer and a scattering length 2«3d. For our ex- 
perimental parameters, we then determine the total number 
of atoms lost due to secondary processes for a thermal sam- 
ple, AOth, and a BEC, Ac^ec- Figure [5] shows that the ra- 
tio (3 + AOth)/(3 + Ao!bec) increases with increasing 030 by 
about 30 percent over the experimentally accessible range of 
fl3D- These results imply that for our experimental parame- 
ters secondary loss processes would result in an increase of 
the ratio g th /g B -^ c in Fig. Ilfb), in contrast to the measured 
data. Thus, within the present model, this rules out secondary 
loss processes as the cause of the effects shown in the present 
work. 



FIG. 5: (color online) Upper estimate for the corrections due to sec- 
ondary loss. For details see text. 



bination event. Here (nl) is the average column density, with / 
the distance covered by a (randomly chosen) atom leaving the 
trap, and a is the scattering cross section. For the atom-atom 



cross section we use the formula a = %na\ D 



/(I 



a„2 



*m) 



where hk is the momentum of the free atom gained in the 
recombination event. For the atom-dimer collision, we use 
a similar expression for the cross section, with the momen- 
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